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1 Introduction 



There has recently been remarkable progress in the understanding of chiral symmetry breaking in 
QCD, based on universality conjectures from random matrix theory g, gj. Central in this develop- 
ment is the idea that the spectra of massless Dirac operators in gauge theories that have non- vanishing 
chiral condensates (as is the case when there is spontaneous chiral symmetry breaking) display uni- 
versal large-volume scaling laws near the origin. There is now mounting evidence that this scenario is 
correct. First, the appropriate microscopic spectral densities, derived from large- ./V matrix ensembles 
with Gaussian weights |], have been shown to consistently reproduce the exact spectral sum rules 
of Leutwyler and Smilga [Ij. Second, the microscopic spectral densities, and in fact all microscopic 
spectral correlators, have been proven to be universal within the given classes of matrix model ensem- 
bles H|. Third, there is now direct numerical support from Monte Carlo simulations of lattice gauge 
theory (for the case of gauge group SU{2) and quenched staggered fermions) that the microscopic 
limit of lattice Dirac spectra is as predicted from random matrix theory || . 

In quenched Monte Carlo simulations it is possible to compute, for finite volume, the lattice eigenvalues 
for massless quarks. In more realistic lattice computations one will keep finite quark masses, and at 
most consider mass rescalings towards smaller values as the volume is increased. So for the practical 
purpose of comparison with lattice gauge theory results, it is imperative that finite quark mass effects 
are under control and understood. From a purely theoretical point of view there is also interest in 
the problem of finite masses. Just as one considers a microscopic rescaling of eigenvalues z it seems 
natural, from the Dirac eigenvalue equation, to consider also the corresponding microscopic scaling 
limit in m jl], M. In view of the above, one might hope that such a double- microscopic limit of the 
Dirac spectrum is universal as well, and computable from random matrix theory. 

The starting point of the following computations are the suggested matrix model ensembles that should 
be relevant for massive Dirac fermions ||, From this we shall derive the pertinent spectral correlators 
(and in particular the spectral densities themselves) for Nf massive fermions in the double-microscopic 
limit where both eigenvalues A and masses m are appropriately rescaled with N. As a straightforward 
by-product of the analysis in ref. ||, we shall prove that all of these double-microscopic spectral 
correlators are universal. We shall also show that the resulting microscopic spectral densities are 
consistent with exact spectral sum rules for massive Dirac operators in QCD and QCD-like theories. 
Such spectral sum rules are normally derived for the case of massless Dirac operators Q, but the 
generalization to massive Dirac operators is not complicated (see also ref. |TJ). 

Our paper is organized as follows. In the next section we introduce the chiral unitary matrix ensemble, 
and use the orthogonal polynomial method to iteratively derive the relevant orthogonal polynomials 
for increasing values of Nf. We prove that in the large- N limit the iVth orthogonal polynomials 

P^ f \\; m\ , . . . , rn? Nj ) for each Nf have universal asymptotic limits near the origin in both A = x/N 2 
and masses m 2 = fi'j/N' 2 . From this it follows that also all microscopic correlators, and the microscopic 
spectral densities themselves, are universal. We also write down a simple set of consistency relations 
for the microscopic spectral densities, which follows from the decoupling of heavy fermions. These 
consistency relations are found to be satisfied by our microscopic spectral densities. In section 3 we 
briefly discuss some of the exact spectral sum rules for massive Dirac operators of QCD-like theories, 
and we verify that our double-microscopic spectral densities are consistent with these exact sum rules, 
which have been derived without the use of random matrix theory. As in the massless case |^, |3|, B], 
this provides strong support to the conjecture that the microscopic spectral densities are universal 
not only within the context of random matrix theory, but indeed are exact expressions also for the 
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full spectral densities of QCD in that limit. Section 4 contains our conclusions, and a proof of the 
general- Nf expressions can be found in the appendix. 



2 The Chiral Unitary Ensemble 

Consider 4-dimensional SU (N c > 3) gauge theories coupled to Nf fermions in the fundamental rep- 
resentation of the gauge group. Assume that a non-vanishing chiral condensate £ = (ipip) has been 
formed. According to the conjectures of ref. M the microscopic spectral density 



C 1' 1 ( ^ 

14— oo V4S V4S 



(1) 



of the Dirac operator can be computed exactly in an ensemble of complex block hermitian (2N+\v\)x 
(2N+\u\) matrices M: 

' ° * ^ . (2) 



M 



W 



The partition function is defined by 



Z 



N f 

dWj[ det (M + im f ) exp 
/=l 



N 



tr V{M 2 



(3) 



Here W is a rectangular complex matrix of size N x (iV + |^|). The integer v is identified with 
topological charge, and the space-time volume V4 is identified with 2iV+|^|, the size of the matrix 
M. The integration measure in eq. (|3|) is the Haar measure of W. From now on consider the sector 
of zero topological charge u, which is the case most relevant for comparison with Monte Carlo data 
of lattice gauge theory. However, the general-^ case can be extracted from the formulas we shall give 
below simply by setting \u\ fermion masses to zero. 



Rewriting the matrix integral (g) in terms of the (positive definite) eigenvalues Aj of the hermitian 
matrix W^W, one gets, after discarding an irrelevant overall factor from the angular integrations p|: 



r oo N ( 

n UA i n(Ai+ 

J0 i=l V / 



m, 



-NV(Xi) 



(4) 



In terms of the standard orthogonal polynomial method, we thus seek polynomials P n 1 (A; m\ , . . . , m Nf ) 
orthogonal with respect to the weight functions 



N f 

^(A) = n(A + m2)e- w ( A ). (5) 



So far we have not specified the potential V(A), which can be parametrized in a quite general way by 

9k 
k 



= E T Xk ■ ( 6 ) 



k>l 

t2< 



As was shown in ref. M, when all rrif = the orthogonal polynomials have, for fixed x = N X 
and t = n/N, a universal limiting behavior. When the polynomials can be normalized according to 
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P^ } (0) = 1 the limit is 



N^oo N 2 



J Nf (u(t)y/x) 

n=Nt '■ J -(u(t)^/2) N f 



with 

The restriction to normalizability according to Pn* f (0) = 1 for n = A^i can be rephrased Q as p(0) 7^ 0, 
so the condition is precisely as expected: the macroscopic spectral density must be non-vanishing at 
the origin. One can then identify 

a = V^l) , p(0) = ^ , (9) 

where a gives the upper limit of the support of p(X). Once universality of the orthogonal polynomials 
has been established, the general proof of universality of all microscopic spectral correlators follows 
as a simple corollary || . 

2.1 One massive flavor 

To generalize these results to the case of massive fermions, we first note that for quenched fermions 
Nf = there is no distinction between the massive and massless cases. Consider next the first non- 
trivial case of Nf = l. We here need polynomials orthogonalized according to 

' d\(X + m 2 ) e- NV ^P t {1) (X;m 2 )P^(X;m 2 ) = . (10) 

Our basic observation is that these polynomials are readily found by expanding in terms of the poly- 
nomials of the quenched case, Pj°\x) by use of Christoffel's theorem Q: 

P^iX; m 2 ) = E^TV(A), (11) 
3=0 h j 

up to an arbitrary normalization condition. By the Christoffel-Darboux formula we can thus choose 
p W (x . m) _ H°\-rn 2 )P^(X) - P^(- m 2 )P?\x) 

Pi (X,m) - XT -2 , (12) 

and fix the normalization subsequently. Note that the r.h.s., despite appearances, indeed is a polyno- 
mial in A. 

Because the new factor (A + m 2 ) in the weight function w(X) can be viewed as a contribution to the 
generic potential V(A) which is subdominant in 1/N, the macroscopic spectral density is unchanged. 
Since p(0) 7^ this means that Pj^ (0;m 2 ) 7^ as well in the large- limit, and by redefining, for 
sufficiently large index i, 

P 4 (1) (A;m 2 ) - P. {1) (A;m 2 )/P 4 (1) (0;m 2 ) , (13) 

the orthogonal polynomials of eq. (|l^) become normalized according to P^\0; m 2 ) = 1. The normal- 
ization constants of course become modified accordingly. 
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In the double-microscopic limit where x = N 2 X and p = Nm are kept fixed as iV^oo we thus get from 
eq. (12) and the previous universal result ([?]): 



" 2 ' Mu(t)VS) + ^J-MWdW'M (14) 



n=Nt 



x + p 2 \ pli{u(t)p) 



where J n (x) and I n (x) are ordinary and modified Bessel functions, respectively. This depends on 
the given potential V(A) only implicitly through the two functions u(t) and r(t), which are inherited 
from the massless case. Because of relation @, we have therefore proved that the double-microscopic 
limit of the polynomial P^\x/N 2 ; p 2 /N 2 ) is universal, depending only on the end-point of the cut a 
(which one can choose at some fixed value), and the value of the macroscopic spectral density at the 
origin, p(0). In applications to Dirac spectra, this value is, by the Banks-Casher relation |10[| ; fixed 



by the chiral condensate: p(0) = (ipip) /V. As a first simple check we note that P^\x/N 2 , p 2 /N 2 ) -► 
J\{u(t)y/x) I '{u(t)yfx) as p— >0, and we hence recover the result of ref. ||. In the opposite limit, p— >oo, 
we see that the above universal polynomial approaches Jo(u(t)y/x), the result of the quenched case. 
This is the first example of the decoupling of a massive fermion, which we shall return to in greater 
generality below. 

It now follows as a simple corollary that also all microscopic spectral correlators are universal in the 
above scaling limit. Consider first the kernel 



Kff( Zl ,z 2 ;m) = yJ\ Zl z 2 \ ^e~ N{ y^) +v ^)) ( Z \ + m 2)( z 2 + m 2) 

X ~2 ~2 ( 15 ) 

z l z 2 

governing the correlation of M . In the scaling limit it becomes of the universal form 

K?((,C';p) = lim 1 A) 
5VS ' s ' p; N->ocN NK N N N J 

f Jo(27rp(0)C) C^i(27rp(0)C) C 2 J 2 (27rp(0)C) 

det Jo(27rp(0)C') (' Ji (2^(0)0 C' 2 J 2 (27rp(0)O 

n( \ KFTii V / o(2vrp(0)^ -p h(2np(0)p) p 2 I 2 (2tt P (0)p) , 
= C(p)^J\(C\ i '- . (16) 

(C 2 -C' 2 )\/(C 2 + » 2 )(C' 2 + » 2 ) 

C{p) is a function of /i only, to be determined below. 
All double-microscopic spectral correlators 



p%\z u ...,z s ;m) = (j] ^tr5(* a -M)\ = det (z a , z b ;m) (17) 



\a=l 



therefore reach universal limits as well: 

p§\Cu • • • , ( s ; p) = Jm • • • , |; £ ) = det K^iCa, C a ; p) . (18) 

In particular, the spectral density itself, 

p$(z;m) = K$\z,z;m) (19) 
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takes a universal form in the double-microscopic limit, which, after choosing the conventional normal- 
ization ||], reads 

pW(C w ) = -L-Kt}\^—,^—;—^—) . (20) 
vs vs>w s v 2 7rp(0)' 27rp(0)' 27rp(0) ; v ; 

After making use of some Bessel function identities, this can, be simplified tcQ 



Ps (C,/*)-y(Jb(C) +Ji(C)J-|CI ( C 2 + ^ )/o(m) 



(21) 



The microscopic spectral correlators, and the density itself, have had their overall normalization fixed 
by the matching condition (the compensating factor of 7rp(0) on the left hand side of the relation 
below is due to the normalization convention (p0|)) between micro- and macroscopic densities, 

limM0)^((;M)] = p(0) , (22) 

the latter of which is independent of p. This is just a convenient short-cut. We could have avoided 
it by explicitly computing the normalization constants hn of eq. ([H]), but which we see no need to 
do that here. In addition we of course have that pg (C'>0) agrees with the result obtained directly 
from the massless case A more interesting condition comes from the decoupling of a very massive 
fermion. This implies 

lirnp«(C; M ) = 4 0) (O- (23) 
One verifies that indeed this relation is satisfied by the density (pip. 

2.2 More massive flavors 

Next, consider the case of more quark flavors. The decoupling of heavy fermions, of which we already 
saw one example in eq. (|23|), leads to a hierarchy of consistency relations which must be satisfied by 
the microscopic spectral densities: 

Ps (C; pi, ■ ■ ■ , PN f ) -> p s (C; mi, • • • , Wf-i) as PN f ^oo 

PS (C; Pi, N f ) -> P S {C;Pl,---,PN f -2> as p Nf -t,p Nf ->• OO 

py\t; Hi, . . . , UNf) -> Ps\0 as pi,...,p Nf -> oo , (24) 

and all other relations obtained by permutations of the pi. There are similar consistency relations 
for all microscopic spectral correlators. In addition, we of course also have that in the limit where all 
masses are set to zero, the result must agree with that of ref. [|J. 

We construct higher- Nf microscopic spectral densities by straightforward iteration of the procedure 
described above. For Nf = 2 the new orthogonal polynomials become 

p(2)a 2 2 . Pj 1] (-mlmDP^X- mf) - iffi (-mj; m\ (A; mf) 

p i '(^,rn 1 ,m 2 ) = , (25) 

A -p //to 



*Our expression for the double- microscopic spectral density does not agree with that of ref. Jj||, where a derivation 
was given in terms of a saddle-point evaluation of a matrix model with Gaussian action. Because of this disagreement, 
we have taken special pains to check that the result presented here is correct. While both the result presented in ref. «n 
and the expression shown here (eq. (|5l])) reduce to the previously known result in the massless limit, only the expression 
reduces correctly to the quenched result in the limit ji — » oo (see below). One further check on the result (^l|) comes 
from the fact that it correctly reproduces an exact spectral sum rule of QCD; see next section. 
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which we again can give our conventional normalization by letting 

P^ 2 \X;mlm 2 2 ) - ff } (A; mf, m|)/^ (2) (0; mf, mi) . 
Taking the large- iV limit with x = N 2 X, n\ = Nm\ and \X2 = Nm2 fixed, leads to 



(26) 



det 



lim pWe^L-A. ifL 

7V™io n { N 2 'N 2 'N 2 ' 



Mi 



m1 



jQ(u{t)y/x) y/x Jl(u(t)y/x) X J 2 (u(t)y/x) 

i (u(t)m) -nxix(u(t)ni) n\h{u(t)m) 



n=Nt 



x + n\ x + ji\ 



det 



-Hih{u{t)ii\) n\h(u{t)n\) 
-[12 h {u(t)fjb 2 ) v\h (■"(*) M2) 



(27) 

Also here, and for the appropriate generalization to all higher values of Nf, the universality proof 

of ref. || immediately implies that Pn^{x/N 2 ; fi 2 /N 2 , y\/N 2 ) is universal. All double-microscopic 
spectral correlators (|l8| ) therefore are universal as well. The results (generalized Bessel kernels) can 
immediately be written down, following the definitions (18). Most compact expressions seem to be 
obtained by writing everything exclusively in terms of the zeroth and first (modified) Bessel functions. 
We restrict ourselves to displaying only the most important of them, the double-microscopic spectral 
density itself, which becomes 



pf&nuto) = ^(Jo(0 2 + Ji(0 2 ) 



(28) 



|C|(/X?-/X|) [/Zl/l(/Xl)Jo(C)+ao(w)Jl(0][M2/l(/X2)Jo(C) + ao(M2)Jl(C)] 



(C 2 +m 2 )(C 2 +m!) mhMhiw) 

We note that it satisfies all the required consistency conditions (E< 



The iterative procedure described above gives the relevant orthogonal polynomials, and hence all 
the microscopic spectral correlators and densities for arbitrary Nf. For example, the orthogonal 
polynomial for Nf = 3 reads 



N^oo n { N 2 'N 2 ' N 2 ' N 2 



(29) 



det 



i4 



n=Nt 
( J (u(t)y/x) 

i (u(t)m) 

Io(u(t)fi 2 ) 
V J («(t)/i 3 ) 



X J\(u{t)y/x) 

ii\h(u(t)m) 



X J 2 (u(t)y/x) 

i4i 2 (u(t)(j, 2 ) 



x 3 / 2 J 3 (u(t)^) \ 

-f4I 3 (u(t)fJH) 

-l4h(u(t)fJ- 2 ) 

-t4i 3 (u(t)n 3 ) ) 



X + fi 2 X + fj, 2 , X + fl 2 




f -mii(u(t)[jLi) 


Hll 2 {u(t)m) 


-nlh (u(t)m) \ 




det 




/J%I 2 (u(t)fJ, 2 ) 


-l4h(u(t)iM 2 ) 






\y -fl 3 Il(u(t)H 3 ) 


f4i 2 ( u (t)fi 3 ) 


-i4h( u {t)ii 3 ) ) 



We shall give the general expression for spectral correlators for an arbitrary number of flavors Nf 

(3) 

below. In fig. 1 we show the microscopic spectral density p s (£; fj,, fj,, y) with, for simplicity, degenerate 
masses. The decoupling of heavy fermions that is expressed by eq. ( p4|) is easily seen in the convergence 
towards the quenched result p% '(C); while in the opposite limit of zero masses the density approaches 
the result of ref. 101 . 
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Figure 1: Microscopic spectral density for 3 flavors with degenerate masses. 



Before proceeding to the general formula for arbitrary Nf, we note that there are still a few particular 
cases that can be expressed in a very simple analytical form. As an example, consider the case of one 
massive fermion (of rescaled mass //), and Nt—1 massless ones. For Nf = 1 this case of course coincides 
with what we have already described above. For general Nf the formula can worked out completely 
analogous to the case of one massive flavor, by starting from the asymptotic form of the polynomial 
(0). The resulting double- microscopic spectral density reads as follows: 



(N f ),, n n ^ \C\j T (n * T ,~ , M 2 4v/-i(C) flNf+M 



One verifies that also this microscopic spectral densities satisfies the consistency condition (pj). 



JiV>-l(0 + ^V>+l(0J 

(30) 



2.3 An arbitrary number of massive flavors 

For practical purposes, only the massive microscopic spectral densities for a few flavors may be of 
relevance. Nevertheless, it is worthwhile to note that the iterative procedure described above can 
be carried through for an arbitrary number of flavors Nf, which we here for notational convenience 
denote by a. 

We shall first state the main results. The orthogonal polynomials, normalized by 

,2 



read 



N 2 ' 



EL 

N 2 ' 



1 



li m p{«)(JL.£. ... ifL 
N^oo n l A 2 'A 2 ' ' N 2 ' 



det 



( J Q {u{t)y/x) 

I Q (u(t)nx) 



n 



,4 



n=Nt 

X Jl(u(t)y/x) 

lixh(u(t)ni) 



\ I (u(t)fi a ) -Hah(u{t)n a ) 



X^' 2 J a+1 (u(t)yffi \ 

{-^) a+1 I a+l {u{t)^) 



{-m) a+1 I a+1 (u(t)fi a ) J 



X + (J, 2 



detij(-fj,iyij(u(t)ni) 



(31) 
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The double-microscopic spectral correlators 



4 a W--,C s ;mr--,/0 = ^^P^ (32) 



are expressed through the kernel which is given by 

ff< a >(C,C; {/*/» = 

/ Jo(2vrp(0)C) C^i(2^p(0)C) 

J (27rp(0)C') C'^l(2^(0)C') 
det Jo(27rp(0)/ii) -^1/1(27^(0)/^ 



(33) 



-vrp(O). 



V^CCI V J (27rp(0K) -^/i(2ttp(0K) 



C a+1 J Q+ i (2ttp(0)C) \ 
r +1 J a+ i(2vrp(0)C / ) 
(-/ii) a + 1 / a+ i(27rp(0)pi) 

(-^) Q+1 / a+ i(2vrp(0)^) / 



C 2 " C' 2 n/=i V / (C 2 + ^/ 2 ) (C' 2 + M/ 2 ) det 4J (-/i 4 )^ 1 / J -i(2vr /3 (0)/i 4 ) 

The double-microscopic spectral density thus explicitly reads 



det 



/ rV-i(C) MO 
MO CMO 
-ft (mi) -Mi h (mi) 



p£°(C;{m/}) 



C%(0 \ 

r +1 J a+ i(o 

(-Ml) Q+1 /a+l(/il) 



ICI V VMa) -llghiPa) ■■■ (-Ha.) a+ I a+ l(Ha) J 

2 n/=i (C 2 + M/ 2 ) detM 



(34) 



Here A4 is an a x a matrix defined by Aiij = (— pj)-' 1 / ? _i(/ij). The proof of the above formulae 
summarized in Appendix A. 



is 



We can now easily check that the microscopic kernel (|33| ) and density (g4|) for arbitrary a satisfies all 
the decoupling relations (p4j) . Suppose p a is taken to be infinity in eq. ([34]) . Due to the asymptotic 
behavior of the modified Bessel functions 



as ijl — > 00 , 



(35) 



the determinant in the numerator of ( p4[) is dominated by its minor for the lower right corner, 



(~Ma) 



a+1 



det 



/ r^-i(o Jb(o 

Jb(C) CMC) 

-Jo(mi) -Mi A (mi) 



\ io(Ma-l) -Ma-l- f l(Ma-l) 

Similarly, its denominator is approximated by 



C^Ja-liO \ 

( a Ja(0 
(-Ml)°/a(Ml) 

(-/i a -l) a Ia(Mai-l) / 



(36) 



0-1 



m. n(c 2 +M/ 2 )-(-M a r- ] 

/=1 



-»/-'a 



det (-^i) J 1 J r j-i(A*i) 



(37) 



\/27rM 

a l5~*)i^a— 1 

Thus we recover the same expression (pij) for a — » a — 1. By iteration, we confirm that the tower of 



consistency relations (24) is satisfied. 
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On the other hand, we have another consistency relation that when all masses vanishes, (£; 0, • • • , 0) 
should agree with the result obtained directly from the massless case ||. In order to check this, let 
us take the limit where the masses are taken to be zero one by one. Suppose fi a is taken to be zero 
first. Since 

/ (0) = 1 , J n (0) = (n > 1) , (38) 

the determinants in the numerator and denominator of ( |34| ) are replaced by their minors for the lower 
left corner. Then we obtain 

/4 a) (C;Mi,---,/«a-i,o) = 



det 



-ICI 



MO 
MO 

-pih(pi) 



\ -p a -ih{p a -i) 



C^-WC) 

cuo 



( a U0 

r +1 J Q+1 (o 



(-(J, a -l) a Ia(}-la-l) {-pa-l) I a +l{Pa-l) j 



2 C 2 U%1 (C 2 + h 2 ) deti< f j< a _i (-myijim) 

We may iterate this procedure as many times as required to obtain, 

( ( a - 2 JaM0 C^JaMO C a U0 

(^JaMO ("MO C +1 J a+ l(0 



(39) 



det 



p^iC^O, •••,0) 



-ICI 



C 2(a-1) ( C 2 +/U 2) (_ M )a-l /Q _ l(M ) 
' C^^a-lCC) ("MO 



(40) 



4 Q) (C;0,---,o) 



-ici det V C%(0 



C Q+1 J a+1 (C) 



c 



2a 



^(j a (C) 2 -J a -i(C)^ + i(0) • 



(41) 

We note that the expression (|40| ) coincides with that of eq. (|30D. The last expression, (|4l]) indeed also 
agrees with the known formula for massless quarks 0. 



3 Exact Massive Spectral Sum Rules in QCD 

So far our discussion has been entirely within the framework of random matrix theory. A highly 
non-trivial question is whether the resulting universal spectral correlators are exact statements about 
the spectral correlators in QCD as well. With the increasing evidence that this is the case for massless 
Dirac operators, it is natural to expect that it extends to massive Dirac operators as well, once 
considered in the appropriate double-microscopic limit which we defined in the preceding section. We 
shall here present some evidence that this is the case. We do this by temporarily leaving the framework 
of random matrix models, and turning to the QCD partition function, as it can be represented in the 
range 

1/Aqcd < N l / d «l/m r Here N gives the size of the volume, and Aqcd is a typical hadronic 
scale in QCD. Following the conventional notation, we denote (up to a sign, see the discussion in 
ref. Q) the chiral condensate by E = ('tpip}- 

Leutwyler and Smilga Q] have shown how to derive exact spectral sum rules for massless Dirac 
operators in the above range|] One generalization to the case of massive spectral sum rules has 

^ See also the work of ref. jllj for generalizations to gauge group SU(2), and to more exotic flavor symmetry breaking 
patterns. 
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briefly been considered by Shuryak and Verbaarschot Since the masses and the eigenvalues are 
rescaled at the same rate in the N — » oo double-microscopic limit, it is clear that we should consider 
spectral sum rules for which the masses and the eigenvalues enter on equal footing. Consider first the 
case of one flavor, Nf = l, in the sector of vanishing topological charge. The natural generalization of 
the spectral sum rules in terms of inverse powers of z 2 is to consider inverse powers of z 2 + m 2 . From 
the expansion 0] of the QCD partition function in terms of Bessel functions in the above range, one 
finds that 

- (42) 



where we have defined p = mNT,. The sum on the right hand side runs over the real zeros jo n of the 
Bessel function Jq{x). Defining the double-microscopic spectral density 

Ps(C^) = ^^(ejv) > p = mNH fixed (43) 

in terms of the ordinary macroscopic spectral density 

p(z) = (£ '8{z - z n )) , (44) 

n 

one sees that the above spectral sum rule can be written 







d( #4 = "Wr • (45) 
C + P 2 2/i/ (/i) V ; 



This massive spectral sum rule, derived entirely within the framework of QCD (a very simple chi- 
ral Lagrangian, in the above range), can now be tested on the double-microscopic spectral density 
p$(Ci p) of eq. (pi]). Performing the required integrals, one finds that it indeed is satisfied. As in the 
massless case, the microscopic spectral density of the chiral unitary ensemble is thus consistent with 
the spectral sum rule of QCD. 

Another example of an exact spectral sum rule which can be derived directly from QCD, and for which 
we can compare with our present results, is given by the case of one massive fermion of mass m, and 
Nf — 1 massless fermions. Again expanding the Bessel functions of the relevant partition function of 
QCD with one massive flavor and Nf — 1 massless flavors [||, one finds the following exact spectral 
sum rule (j m , n denotes the nth zero of the Bessel function J m (x))\ 



C 2 + V 2 n 3Nf-l,n + V 2 

2p,I Nf -i{n) 
, . . . , 0) of eq 

is satisfied. Other spectral sum rules can be worked out analogously. 



(46) 

tpj-Nf-lKP) 

Inserting the microscopic spectral density Pg {& p,0, . . . ,0) of eq. (|3^) we verify that also this identity 



10 



4 Conclusions 



Motivated by the recent progress in numerically computing the microscopic spectral densities of the 
Dirac operator in realistic four-dimensional lattice gauge theories ||, we have set out to compute 
the microscopic spectral densities of what should correspond to massive Dirac operators. Such an 
extension will presumably be required before a detailed comparison with lattice gauge theory data 
can be performed beyond the quenched approximation. To reach the universal limit, masses must be 
scaled at a very precise rate as the volume is increased, but this is entirely feasible in the context of 
lattice gauge theory. We thus expect that the results presented here may be of practical value when 
it comes to detailed comparisons between lattice Monte Carlo data for QCD and the universality 
predictions from random matrix theory. 

We have succeeded in deriving the relevant microscopic (we have called these "double-microscopic" 
because both masses and eigenvalues are rescaled as N — > oo) spectral correlators, and in particular 
the microscopic spectral densities themselves, within the framework of random matrix models. In 
doing so, we have simultaneously extended the universality proofs of ref. || to this more general 
situation.^] The case of the unitary matrix ensemble, conjectured to be possibly relevant for SU(N C ) 
gauge theories in (2+1) dimensions ||, can be worked out analogously and will be presented elsewhere 
|14j |. As in ref. ||, the proven universality is strictly limited to the framework of random matrix 
theory. 

Based on the decoupling of heavy fermions, we have derived a set of consistency conditions for the 
microscopic spectral densities. These are non-trivial relations that show how the various microscopic 
spectral densities Pg* f \C', Mi> • • • , A 4 at* ) for different values of Nf must be related to each other as one 
or more of the masses are sent to infinity. We have verified that our spectral densities satisfy all of 
these general consistency conditions. In the other extreme limit where all masses are taken to zero, 
we of course recover the known results [||, [|. 

To investigate the question as to whether QCD with massive flavors indeed fall into the universality 
classes derived here, we have confirmed that exact massive spectral sum rules (generalizations of 
the Leutwyler-Smilga sum rules Q |l|) derived directly from QCD are satisfied if we identify the 
double-microscopic spectral density of QCD with that of the chiral unitary ensemble. This, together 
with the mounting evidence from the massless case, makes it highly plausible that all the universal 
double-microscopic spectral correlators we have derived here are exact statements about QCD in the 
above limit. The crucial test will be a comparison with results from lattice gauge theory, where these 
microscopic spectral correlators and densities are very convenient finite-size scaling functions, almost 
tailored for Monte Carlo simulations. Hopefully such data will soon be available. 
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Note added: Our main results (33) and (34) of this paper have been derived independently for 

*An alternative route would be to derive directly the differential equation that the orthogonal polynomials satisfy in 
the microscopic limit, for example using the method of Kanzieper and Freilikher p^| , or to use the addition formalism 
ofZee(l|. 
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the case of the Laguerre ensemble by Wilke, Guhr and Wettig ( hep-th/9711057 ), in a paper which 
appeared on the e-print archives just a few days after ours (these authors did not consider the question 
of universality). The sum rule (45) has been derived earlier by Gasser and Leutwyler in Phys. Lett. 
B188 (1987) 477, and has been compared to lattice QCD data by Verbaarschot in Phys. Lett. B368 
(1996) 137. The authors of ref. j7| have found an error in their previous calculation of the Nf = 1 
case in the Gaussian ensemble, and they now agree with the expression given in eq. (23) (M. Nowak, 
private communication) . 



A Proof of the Formulae in section 2.3 

Since the procedure of constructing the kernel out of the polynomials (|l5|) is essentially the same as 



that of adding another massive flavor (12) up to the sign of m 2 , it suffices to prove (33). We first 
demonstrate the following lemma: 

Lemma : 

Let P^ a \t] Ao, Ai, • • • , X a ), a = 0, 1, 2, • • •, be a set of functions generated by the iteration 

p[ Q+1 ](t;A ,Ai,---,A Q+1 ) = (A.l) 
P^(t; Ao, Ai, • • • , \ a )pj; a \t; X a+ i,X 1 • • • , A a ) — P t a (t; Aq, Ajv • , X a )P^(t; A a +i, Ai • • • , X a ) 



Aq — A, 



Q + l 



Then they are given by 



P*- a '(t; Aq, Ai, • • • , Aq,_i, A q ) — c(t; Ai, • • • , X a -i 



detij P^ (t; Xj) 



(A.2) 



where p(*)(t; A) = ^P [0] (t; A), and c(t; A x , • • • , A a _i) is a function in t and in {Ai, • • • , X a -i}. 



Proof : We prove it by induction. Suppose ( |A.2j ) holds for an a. Let PW(A) = dLpl°}(p,X) (the 
argument t is for convenience suppressed below). Consider the numerator in (A.l), 



p p> ... p(a-l) p(«)(Ao) 



c(Ai, • • • , A a _i) 

— (Ao «-> A a+ i) 

= c(Ai, • • • , A Q _i) 2 

— (Aq «-> A a+ i) 



P P' 



p(o-l) pW(A 1 ) 



/ 



p p> ... p(a-l) p( Q )(A Q ) 



P p> ... p(«-i) p(°)(A ) 
p p> ... p(a-l) p(«)(Ai) 

P p ... p(«-i) p(«)(A a ) 



c(Ai, • • • , A Q _i) 



V 



p P' 
p P' 



p p ... p(a-l) p(«)(A a+1 ) 
p p ... p(a-l) F(«)(Ai) 

p p ... p(a-l) P( Q )(A Q ) 



p(a-l) p(a+l)(A a+1 ) 
p(a-l) p( Q+1 )(Ai) 



p p ... p(«-l) p( a +!)(A a ) 



c(Ai, • • • , A a _i) 2 P(Aq, Ai, • • • , A a _i, Aq) . 



(A.3) 



The above identity holds because the t— derivatives acting on c(Ai, • • • , A a _i) produce terms indepen- 
dent of Aq and A Q +i and are thus cancelled by the subtraction interchanging Aq «-> Aq,+i, and those 
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acting on the columns 

p(a-l) 

vanish by themselves (the same columns are already present in the determinant). 

By definition, D(Xq, Ai, • • • , A a _i, A a ) enjoys the following properties: 

• It is a sum of monomials of the form 

± P(#)(A )P( # )(A 1 ) J P(#)(A 1 )...P(#)(A Q )P(#)(A Q )P(#)(A Q+1 ) , 
with (total # of derivatives) = 2(0 + H \- a) + 1 = a 2 + a + I. 

• It vanishes if any two A's are coincident. 

• It is antisymmetric under Ao \ a +i- 

• It is completely symmetric under Aj <-> Xj, 1 < i, j < a. 

The only function which satisfies these properties is, up to a constant, 







p 


P' •• 


. p(<*- 


-i) 


p{a) 


P (a+1) (Ao) 




P p' ... p( Q - 1 )(A 1 ) 




p 


P' •• 


. p(a- 


-i) 


p(a) 


p(«+i)( Al ) 




P p' ... p( a -!)(A Q ) 




p 


P' •• 


. p(a- 


-i) 


p(a) 


P (Q+1) (Aa) 






p 


P' •• 


. p(a- 


-i) 


p(a) 


P(« +1 )(A Q+1 ) 


On the other hand, the denominator in 


(A.l) 


is given 


by 








(A -A 


a+i) n (Aj — n 


(A; 


- A, ) = 


n 


(\ 


-A,) 


n 



0<i J <a 



0<i,j<a+l 



(A.4) 



(A.5) 



(A.6) 



Since the first factors in (A~5) and (A.6) as well as c(Ai, • • • , A a _i) 2 are independent of Ao and X a+ ±, 
it can be absorbed into the redefinition of c(Ai, • • • , X a ). Then we recover (|A.2 ) for a + 1. QED. 

Now we replace A, -> ( 2 /N 2 and pM(t,A) by its microscopic limit, P^(t; ( 2 /N 2 ) -> J (u(t)(). Then 
we can inductively prove that its ^derivatives are expressed as 



(A.7) 



k=0 



Once it is substituted inside the determinant det P^ (Xj), only the top term proportional to C Ji(u(t)C) 
contributes. Thus the determinant in (A. 2) is replaced by 



IIM*) det CJi(n(t)C 



-i=0 



0<«,i<a+l 



(Ai 
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In order to construct the kernel, we make an analytical continuation of (Ci>'"'>Ca) to imaginary 
(ifj-l, • • • , i^ a ), denote the rest as Co = C> Ca+i = C'> an d multiply the whole expression at t = 1 

by the factors from the measure and the Jacobian, J\(('\ ]l/=i(C 2 + M/XC' 2 + M/)- Then we obtain 



^)(C,C / ;{M) = ^(W})v / ICC 





Jo(27rp(0)C) 


CJi(2vrp(0)C) ■■ 


C a 


^J a+ i(2vrp(0)C) 




Jo(2ttp(0)C') 


C Ji(27rp(0)C) •• 




fl J Q+1 (2vrp(0)C / ) 




J {2tv P (0)m) 


-/x 1 / 1 (27rp(0)/xi) •• 




Q+1 / Q+ i(27rp(0)/zi) 


'icci- 


I (27rp(0)yU a ) 


-lu, Q Ii(2irp(0)n a ) ■■ 


• (-Ma) 


Q+1 / a+ i(27rp(0)^) 



(A.9) 



/4 Q) (C;{M) = C({M)KI 



C-^-i(C) 


MO 


• C a -^a-l(C) 


C a UC) 


Jb(0 


CMC) ■■ 


C%(0 


C a+1 J«MC) 


^o(Mi) 


-m Mm) •• 


• {-mTUm) 


(-M) a+1 iaMM) 


Jo (Ma) 


-HaMHoi) 




(-fi a ) a+1 i a+1 (fi a ) 



n/=i(c 2 + ^) 



(A.10) 



The constant C({n f }) = (27rp(0)) (1 "t)( 1+a ) C({/x/}) is determine d to be -(2detA4) _1 by requiring 
the matching between the £ — ► oo limit of the microscopic density ( A.lCj ) and the macroscopic density: 



det M 



p ( sHc 



oo 



{H f }) - C({ M/ })|C|- 



( a - l J a M() CMC) 
C a J a (C) C a+1 J*MC) 



c 



2a 



-C({n f })-detM 

TT 



(A.ll) 



By our normalization (|2Q), this should equal l/ir. Substituting the resulting C back to ( |A.S| ) and 
( |A.1C| ), we establish the result announced in section 2.3. 
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